The Phase Structure of the Polyakov-Quark-Meson Model 
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The relation between the deconfinement and chiral phase transition is explored in the framework 
of an Polyakov-loop-extended two-flavor quark-meson (PQM) model. In this model the Polyakov 
loop dynamics is represented by a background temporal gauge field which also couples to the quarks. 
As a novelty an explicit quark chemical potential and iV/-dependence in the Polyakov loop potential 
is proposed by using renormalization group arguments. The behavior of the Polyakov loop as well as 
the chiral condensate as function of temperature and quark chemical potential is obtained by mini- 
mizing the grand canonical thermodynamic potential of the system. The effect of the Polyakov loop 
dynamics on the chiral phase diagram and on several thermodynamic bulk quantities is presented. 
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I. INTRODUCTION 

Driven by the heavy-ion programs at GSI, CERN SPS, 
RHIC and soon the LHC there is strong interest in the 
properties of strongly interacting matter at extreme tem- 
peratures and baryon densities. Ultimately these have to 
be understood on the basis of Quantum Chromo dynam- 
ics (QCD), which governs the strong interaction sector of 
the Standard Model. 

QCD at zero temperature and density is well- 
established by now both numerically and analytically. 
Besides numerical evaluations on discrete space-time lat- 
tices, e.g. [H, 01i functional methods based on the Func- 
tional Renormalization Group (FRG), e.g. 0, Q], and 
Dyson-Schwingcr equations (DSE), e.g. [5|, have been 
used to elucidate our understanding of the theory of 
strong interactions. If applicable, lattice computations 
give high numerical accuracy with small truncation er- 
rors. Functional methods have their advantages if it 
comes to the deep infrared, the simple explanation of 
physical mechanisms as well as the inclusion of chiral 
dynamical quarks. The last years have seen a very fruit- 
ful interaction between the different methods leading to 
a largely quantitative understanding of QCD at vanish- 
ing temperature and density even though the full undcr- 
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standing of the confinement mechanism and its relation 
to spontaneous chiral symmetry breaking is not yet set- 
tled. 

At finite temperature and in particular at finite baryon 
density or chemical potential fig the situation is much 
less clear. At finite /is lattice computations have princi- 
pal limitations due to the complex action which hampers 
stringent theoretical evaluations of the QCD phase dia- 
gram, in particular the possible critical endpoint of the 
line of first-order transitions from first principles. Func- 
tional methods have been used to obtain results for pure 
Yang-Mills at finite temperature, as well as the hadronic 
sector of QCD at finite temperature and density, see 
e.g. [1, 0|- A full QCD study is hampered by the fact 
that the gauge sector, i.e. the confinement-dcconfincmcnt 
phase transition in pure Yang-Mills, is not fully resolved 
yet: the potential for the order parameter, the Polyakov 
loop, does not lead to a phase transition in perturbative 
computations, e.g. but recently this gap has been 

closed within a non-perturbative flow study [9(. Hence, 
a first principle approach to QCD at finite temperature 
and density with functional methods is in reach. In our 
opinion this opens a way to continuing the fruitful in- 
teraction between the different methods that has already 
proven so successful at vanishing temperature. 

A first step towards a full QCD study with functional 
methods is done by studying effective Lagrangian models 
which are constructed from the non-perturbative Yang- 
Mills effective potential and effective hadronic models. 
In recent years a promising realization of this idea has 



2 



been put forward, based on lattice results for the ther- 
modynamic potential in pure Y ang -Mil ls ( YM) theory 
and universality arguments [lol . I 111 . E2I [l3l . 
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Pure YM theory corresponds to the 
heavy-quark limit of QCD in which the Polyakov loop ex- 
pectation value serves as an order parameter for confine- 
ment. This approach results in an effective scalar Z^ c - 
theory whose physical minima are the vacuum expecta- 
tion values of the Polyakov loop. It has been observed 
that a thermodynamic potential for the Polyakov loop 
can be constructed where the parameters are fitted to 
precise finite-temperature lattice data for the equation of 
state (EoS) in the heavy-quark limit. This very successful 
theory simulates the calculated first-order confinement- 
dcconfincmcnt phase transition at finite temperature. 
For two light flavors, on the other hand, QCD exhibits 
an (almost) exact chiral symmetry and is believed to be 
in same universality class as the 0(4) model 22| . An 
effective realization of this symmetry is provided by the 
Nambu-Jona-Lasinio (NJL) model or, upon bosonization, 
the quark-meson (QM) model. It is therefore natural to 
combine both aspects of QCD by coupling light chiral 
quarks to the Polyakov loop field. This results in the 
PNJL model 0, Q El 0, 0] or the PQM model, which 



has the benefit of renormalizability, and a simpler linkage 
to full QCD, see e.g. 25|. Calculations for the thermody- 
namic potential and the resulting phase structure will be 
performed at the mean-field level as was done in similar 
analyses using the PNJL model. Eventually, however, we 
plan to include fluctuations using RG-techniques [3]. 

Even at the mean-field level there remain some open 
issues. First of all, the Polyakov loop potentials sug- 
gested so far, are fixed at vanishing [Lb- In this case, the 
expectation value <J> of the Polyakov loop operator and 
that of its adjoint, $ are linked by complex conjugation. 
At finite chemical potential this relation is lost and the 
effective potential depends on the two independent vari- 
ables $ and Extensions to finite /j,b have therefore 
to be carefully evaluated. Moreover, the flavor and den- 
sity dependence of the pure Yang-Mills potential has not 
been explored as yet. We will show that it is possible 
to extract the qualitative behavior in the fully coupled 
system by perturbative arguments, as well as physical 
consistency arguments. 

The outline of the paper is as follows: in the next 
section wc introduce the Polyakov loop variable and dis- 
cuss its effective potential derived from lattice data. The 
Polyakov loop potential is then coupled to the quark- 



meson model which defines the Polyakov-quark-mcson 
model. The grand canonical thermodynamic potential of 
this model is derived in mean-field approximation, and 
the choice of model parameters is discussed. Sec. IIIII 
is devoted to thermodynamical applications, in partic- 
ular we evaluate the pressure, quark number density and 
quark number susceptibility. Moreover, the phase struc- 
ture, i.e. the chiral and confinement-dcconfincmcnt phase 
transition is explored. Subsequently, the influence of the 
Polyakov loop potential on the thermodynamics is inves- 
tigated and in Sec. IIVI concluding remarks are drawn. 



II. POLYAKOV-QUARK-MESON MODEL 
A. Polyakov loop potential 

A key observable in QCD at finite temperature is the 
Polyakov loop. Its expectation value serves as an order 
parameter for confinement in the heavy-quark limit. The 
Polyakov loop operator is a Wilson loop in the temporal 
direction and reads 



V{x) = "Pexp 



i j cItAo(x,t) 



(2.1) 



where V stands for path-ordering and Aq (x, t) is the tem- 
poral component of the Euclidean gauge field 
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30| | . The color trace of (|2.ip in the fundamental rep- 



resentation tr c "P(:zf) is the creation operator of a static 
quark at spatial position x. Periodic boundary condi- 
tions ensure gauge invariance of Eq. (|2.4[) up to center 
elements. This goes hand in hand with the fact that 
the temporal or Weyl gauge Aq = cannot be achieved 
for periodic boundary conditions. Its physical interpre- 
tation is best seen in Polyakov gauge, where the tem- 
poral component of the gauge field is time-independent, 
A(\{x, t) = ^4 q (x) , and is in the Cartan sub-algebra (see 
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34[). Hence, within this gauge the 



eg- 

Polyakov loop operator simplifies to 
V(x) = cxp (i/3A c {x)) 



(2.2) 



with A c {x) = A^\x)t 3 + ^ 8) (x)r 8 . This results in a 
simple relation between the Polyakov loop and the tem- 
poral component of the gauge field, 



Al{S) = -i{d p V{£))V\S). 



(2.3) 



The normalized Polyakov loop variable <j>(x) and its 
hcrmitian (charge) conjugate <£>(af) are defined as the 
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thermal expectation value of the color trace of the 
Polyakov loop operator (|2.1[) 

*(a?) - ±-(tr c V(£))/3 , *(£) = ^-(tr c pt(f) )/3 . (2.4) 

We emphasize again that the traces are taken in the fun- 
damental representation. <&(:?), $(x) are complex scalar 
fields. Their mean values, i.e. the solution of the quantum 
equations of motion, are related to the free energy of a 
static, infinitely heavy test quark (antiquark) at spatial 
position x. The order parameter <S>(x) vanishes in the 
confined phase where the free energy of a single heavy 
quark diverges. In the deconfincd phase it takes a finite 
value. The correlation function of two Polyakov loop 
variables is related to the free energy F q q of two color 
sources q and q with spatial separation r '= x — y as 
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— {tr c V{x)tr c V\y))p = e 



(2.5) 



The dependence on r allows one to extract the string 
tension. The cluster decomposition property (locality) 
enforces that for infinite distance the correlation between 
a quark and anti-quark vanishes and we arrive at 
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(tr c V(x)tr c V\y))p -> *(z)*(y) 



(2.6) 



These properties provide the Polyakov criterion of con- 
finement at finite temperature. It is linked to the cen- 
ter Z/v c symmetry of the SU(N C ) gauge group: a gauge 
transformation that is periodic up to a center element, 
leads to 



$ — > z$, z G Z 



(2.7) 



Thus, the confining phase is center symmetric, whereas in 
the deconfincd phase center symmetry is spontaneously 
broken. 

In summary, the confincmcnt-dcconfinemcnt phase 
transition is characterized by the mean value $ = in the 
confined phase and a finite non-zero value in the decon- 
fined phase. In the presence of dynamical quarks, the free 
energy of a quark-antiquark pair does not diverge any- 
more, and the order parameter is always non- vanishing. 
For finite quark chemical potential the free energies of 
quarks and antiquarks are different. Since $ is related 
to the free energy of quarks and the hermitian (charge) 
conjugate $ to that of antiquarks, their modulus in gen- 
eral differs, i.e. $ ^ In pure Yang-Mills theory the 
mean values $, * are given by the minima of the effective 
Polyakov loop potential U(Q, *). It can be constructed 



from lattice data for the expectation values <f>, $ fl"o| . 
Here we use a polynomial expansion in $, $ up to quar- 
tic terms. This leads to an effective potential U in terms 
of the moduli |$| and |l>|, the product $$, and in $ 3 , <l 3 
related to the Z% symmetry. The £/(l)-symmetric part 
of Li is a Ginzburg-Landau type potential. 

In pure Yang-Mills we can restrict ourselves to fields 
with the same modulus, |$| = |$|. With this additional 
constraint we have = |$| 2 and we can drop one of the 
[/(l)-invariants in the expansion. This has been used in 
Ref. where the potential is only expanded in $*, 
$ 3 and $ 3 . Alternatively one can use the moduli and 
drop the $$-tcrm. We conclude that the general effective 
potential in this approximation reads as 
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1*1 



1*1 



-|($ 3 + $ 3 ) + |(|$| 



1*1 



(2.8) 



The expansion coefficients are fixed to reproduce ther- 
modynamic lattice results for the pure YM sector as in 
Refs. 0,0,SHE3,Eil- This leads to temperature- 
independent coefficients 63 = 0.75 and 64 = 7.5, and a 
temperature-dependent one 62 with 



( y 



'T 



where ao = 6.75, a\ = —1.95, a 2 = 2.625, 03 = —7.44, 
63 = 0.75 and 64 = 7.5. The effective potential (|2 . 8[) can 
be augmented by logarithmic terms, see e.g. Ref. [3^]. 
This will be discussed in Sec. Ill Dl 

The potential (|2.8p with the above parameters has a 
first-order phase transition at the critical temperature 
T = 270 MeV. 



B. Coupling to the quark-meson sector 

The hadronic properties of low-energy QCD with light 
flavors are effectively incorporated by a chiral quark- 
meson model. Here the local SU(N C ) gauge invariance of 
the underlying QCD is replaced by a global symmetry in 
the original quark-meson model which results in the loss 
of the confinement property. The QM model shows a chi- 
ral phase transition at realistic temperatures e.g. [3. [111]. 
In the limit of massless quarks the order parameter of the 
chiral phase transition is the quark condensate (qq) . For 
realistic up- and down quark masses chiral symmetry is 
broken spontaneously and also explicitly in the vacuum 
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resulting in a finite chiral condensate (qq). Due to the 
lack of confinement in this model single quark states are 
already excited at low temperatures in the chirally bro- 
ken phase, see e.g. 40[ resulting in an unrealistic EoS 
near the phase transition. Since the constituent quark 
masses are much larger than that of the pion the meson 
dynamics dominates at low temperatures and the predic- 
tions from chiral perturbation theory are reproduced. 

By combining the Polyakov loop model with the QM 
model chiral as well as confining properties of QCD are 
included. This promi sing appro ach has been put forward 
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proved the EoS near the phase boundary. The integra- 
tion over the gluonic degrees of freedom in the presence of 
a homogeneous background for the temporal component 
Aq yields the Polyakov loop potential and the mesonic 
terms of the chiral QM model. Thus, the dynamical 
quark sector of QCD is included by integrating out the 
quarks in the presence of mean background fields. This 
finally leads to a coupled Polyakov-quark-meson model 
with an interaction potential between quarks, mesons 
and the Polyakov loop variables $, To leading loop 
order this potential is provided by the Dirac determinant 
in the presence of the mean fields. 

The generalized Lagrangian of the linear QM model 
for Nf = 2 light quarks q = (u, d) and N c — 3 color de- 
grees of freedom coupled to a spatially constant temporal 
background gauge field reads 

£ = q(ip-g(a + i l5 TTT))q+^{d^) 2 + ^(d^) 2 

-U(a,Tr) , (2.10) 

where the purely mesonic potential is defined as 

[/(cr,7?) = ^(cr 2 + TT 2 - V 2 ) 2 - «7 . (2.11) 

The isoscalar-scalar a field and the three isovector- 
pseudoscalar pion fields 7? together form a chiral vector 
field <j). Without the explicit symmetry breaking term c in 
the mesonic potential the Lagrangian is invariant under 
global chiral £[7(2) j, x SU(2)r rotations. The covariant 
Dirac operator = — iA^ in (|2.10p reads 



(2.12) 



C. Polyakov loop potential parameters 

In the presence of dynamical quarks, the running cou- 
pling a is changed due to fermionic contributions. In 
our approximation to the Polyakov loop potential this 
only leads to a modification of the expansion coefficient 
&2, Eq. (|2.9j) . The size of this effect can be estimated 

H 
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within perturbation theory, see e.g. 
At zero temperature it leads to an ^-dependent decrease 
of Aqcei > which translates into an Ay-dependent decrease 
of the critical temperature T at finite temperature. The 
two-loop /3-function of QCD with massless quarks is given 
by 



/3(a) = -ba 2 
with the coefficients 

b = i-(HJV c -2JV», 



car 



(2.13) 



(2.14) 



24tt 2 



34 A 2 - lON c Nf - 3- 



N 2 



Nr. 



N f . (2.15) 



Here, we have assumed a RG scheme that minimizes 
(part of) the higher-order effects. This is an appropri- 
ate scheme for our mean-field analysis. At leading order 
the corresponding gauge coupling is given by 

no 



a(p) 



OK) 



(2.16) 



1 + a b\n(p/A) 

with ao = ct(A) at some UV-scale A, and Aqcd = 
Aexp(— l/(ao&)). At p = Aqcd the coupling (|2.16|) ex- 
hibits a Landau pole. At finite temperature the relation 
(|2.16[) allows us to determine the Aj-dependence of the 
critical temperature To {Nf). For Nf = it is given by 
To = 270 MeV which corresponds to fixing the coupling 
ao at the r-scale T T = 1.770 GeV and a running coupling 
of ao = 0.304 accordingly. If one keeps the coupling ceo 
at T T fixed, this identification yields the relation 



T (Nf) = T T e- 1 '^ 



(2.17) 



and Table U for the A/-depcndcnt critical temperature 
Tq in the Polyakov loop potential for massless flavors: 



N f 





1 


2 


2 + 1 


3 


T [MeV] 


270 


240 


208 


187 


178 



TABLE I: Critical Polyakov loop temperature To for Nf mass- 
less flavors. 



The spatial components of the gauge fields have vanishing 
background i.e. A^ = S^qAq. 



Massive flavors lead to suppression factors of the order 
Tq/(Tq + m 2 ) in the /3-function. For 2 + 1 flavors and 
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a current strange quark mass m s w 150 McV we obtain 
T (2 + 1) = 187 McV. We remark that the estimates for 
T (Nf) have an uncertainty at least of the order ±30 
MeV. This uncertainty comes from the perturbative one- 
loop nature of the estimate and the poor accounting for 
the temperature effects. For example, with the two loop 
coefficient (|2.15j) and concentrating on Nf = 2 as studied 
in the present work we are led to To (2) = 192 MeV. 
Fortunately, the results only show a mild Tq dependence. 

Finally, we argue that there are no double counting ef- 
fects due to the inclusion of the Dirac determinant in the 
PQM and the independent adjustment of the Polyakov 
loop model parameters: the Polyakov loop potential pa- 
rameters, in particular 62, Eq. (|2.9|) . genuinely depend on 
the running coupling, which is changed in the presence of 
quarks. This effect is modeled by changing To — > To(Nf) 
as defined in Eq. (|2.17|) . The direct contributions to 
the Polyakov loop potential which originate from the 
fcrmionic determinant f2q q , Eq. (|3.24p , are not governed 
by this redefinition, and have to be added separately. 

D. Non- vanishing chemical potential 



A further intricacy concerns the Polyakov loop poten- 
tial at finite chemical potential 5l|, 152J . Then the con- 
straint i> = <f>^ ceases to be valid, and the extension of 
Eq. (|2.8p to finite /x is not unique anymore. For further 
details see e.g. Refs. 



53 



54j . The leading /x-dependence 
of the full potential stems from the Dirac determinant, 
and we assume that the $, ^-symmetric form of the po- 
tential (|2.8|) persists at finite \x. Then the only additional 
/x dependence originates from a possible /x dependence of 
the model parameters. This approximation certainly is 
valid for small chemical potential where the /x dependence 
is rather small. The remaining ambiguity concerns possi- 
ble <&3>-terms, that can be incorporated into the potential 



2.8[) by the replacement 



m 2 + m 2 ) 



|$| 2 + |$| 2 ) + (1 - 0)$<f>. (2.18) 



Eq. (|2. 18|) leaves the potential unchanged for /1 = 0. 
that is § = <&t, p or positive 8 the potential has un- 
stable directions, e.g. for vanishing $ or $, and large <&, 
$ respectively. Hence the choice 8 = possibly leads 
to negative susceptibilities. In Refs. 18|, [20( this choice 



has been used, and the computed susceptibilities are not 
This problem has been cured in 

Refs. ' 
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39j,|5l| by augmenting the Polyakov loop poten- 



changing the model parameters in the polynomial ansatz 
used in these works. For 8^0 these logarithmic terms 
are not necessary, due to lack of unstable directions. Fur- 
thermore, a weak total /x dependence as well as the valid- 
ity of the mean-field analysis would hint at the preferred 
choice = 1. However, for this choice the expectation 
value of ($— <3?) has the wrong sign, even though other ob- 
servables show a mild 0-dependence. Clearly, this struc- 
ture is related to the present mean-field approximation. 
It should be possible to overcome this parameter depen- 
dence in a fully non-perturbative setting. Here, we shall 
show results for the choice 8 = 0. 

In a final step we implement a /x-dependent running 
coupling in the &2 coefficient, analogous to the Nf- 
dependence discussed above. Indeed, one can argue 
that this is a minimal necessary generalization: without 
a /x-dependent bi the confinement-deconfinement phase- 
transition has a higher critical temperature than the chi- 
ral phase transition at vanishing chemical potential. This 
is an unphysical scenario because QCD with dynamical 
massless quarks in the chirally restored phase cannot be 
confining since the string breaking scale would be zero. 

As for the A^-dependence we resort to perturbative es- 
timates. To begin with we simply allow for an additional 
xt-dependent term in the one-loop coefficient b, 

b(j i ) = ±(UN e -2N f )-b li £ s . (2.19) 

This specific simple choice of the /x-dependent part can 
be motivated by using HDL /HTL results on the effective 
charge [55 1 



a(p, T, /x) = 



a(p) 



1 + m 2 D /p 2 ' 



(2.20) 



tial with logarithmic terms. Effectively, this amounts to 



with the perturbative Debye mass m 2 D = (N c /3 + 
N f /6)g 2 T 2 + N f /(2Tr 2 )g 2 [i 2 . The /x-derivative of the 
modified coupling, [id^a = b^f^/p 2 , can be related to 
a momentum derivative pd p a = —b(p,fi)a 2 . Within 
the present simple approach based on a /x-dependence 
only valid in the perturbative regime we estimate 
the momentum-dependent coefficient b(p, /x) by 6(/x) = 
b{p/T T ,fjb) at an (average) momentum scale 7TV with 
7 < 1. 

The coefficient b^ can be fixed such that the chiral tran- 
sition temperature and the confinement-deconfinement 
transition agree at some non-vanishing /x. Interestingly, 
it turns out that then the transition temperatures agree 
for all /x's. The related value of 6 M is provided by 7 ~ 1/4 
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and 
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-N f 



(2.21) 



Inserting the //-dependent coefficient b(/i) into Eq. (|2.17|) 
then leads to an additional //-dependent To, 



T (ji,N f ) =T T e- 1 l { - olob ^ . 



(2.22) 



Eq. (|2.22p with (|2.2ip should be viewed as a rough esti- 
mate of the //-dependence of To. We emphasise again that 
this simple estimate leads to coinciding phase boundary 
lines for the chiral and confinement-dcconfinemcnt tran- 
sition, see Sec. IIIII For more quantitative results the 
non-perturbative running of the coupling in the presence 
of finite temperature and quark density has to be consid- 
ered. This can be incorporated in a self-consistent RG- 
setting. Moreover, one has to resolve the uncertainties, 
discussed at the beginning of this section, concerning the 
form of the effective potential at finite //. 

Here we will present a comparison of the phase diagram 
with and without //-dependent To in Fig. [6j For the other 
results the additional //-dependence is taken into account. 

III. APPLICATIONS 



The PQM model is defined by the Lagrangian (|2.10D 
with the Polyakov loop potential Eq. (|2.8|) . The depen- 
dence of the coefficient 62 Eq. (|2.9[) on the 7V/-dependent 
or (Nf , /^-dependent running coupling a is governed 
by Eq. ([2~T7]) and ([2~2"2"]) respectively. This defines the 
starting point for an investigation of the phase struc- 
ture and bulk thermodynamics of the PQM model. The 
thermodynamics is characterized by the grand canonical 
potential which is analyzed in mean-field approximation. 

A. Grand canonical potential 

The grand canonical potential in a spatially uniform 
system is determined as the logarithm of the partition 
function, which in our case is a path-integral over the 
meson and quark/antiquark fields including the Polyakov 
loop. We confine ourselves to the SU (2) /-symmetric case 
and set /1 = fi u = fi^. This is a good approximation to 
the realistic case since flavor mixing in the vector chan- 
nel is small. Integrating over the fermions by using the 
Nambu-Gor'kov formalism and introducing averaged me- 
son fields yields the grand canonical potential 



U(^,^) + U(a)+n qq (^,^,a) (3.23) 



with the quark/antiquark contribution 

ln(l + 'P t e- (i? '' + ' i)/T )} (3.24) 

and the purely mesonic potential 
A 



U(a) 



2\2 



ca 



(3.25) 



The divergent vacuum part in the quark/anti-quark 
contribution is absorbed in the renormalization which 
is done in the vacuum. The quark/antiquark singlc- 
quasiparticle energy is given by 



P 2 



(3.26) 



with the constituent quark mass m q = ga. The remain- 
ing color trace in the quark/antiquark contribution (|3. 241) 
is evaluated by using the identity Tr In A = In det A and 
yields 



(3.27) 



{ln[l+3(* + ^e-^-^' T )e-^-^l T +e-^ E --^' T 

l + 3($ + $ e -( £ P+A')/T) e -(£;p+M)/T +e -3(£; p +M)/T 



In 



} 



Note, that no ultraviolet cutoff is necessary because the 
PQM model is renormalizable in contrast to the PNJL 
model (see e.g. [18l. Il9j). 

The equations of motion are obtained by minimizing 
the thermodynamic potential (|3.23[) w.r.t. the three con- 
stant mean fields a, $ and 



0. 



(3.28) 



= ( CT ), *=<*>,*=(*) 



The solutions of these coupled equations determine the 
behavior of the chiral order parameter (a) and the 
Polyakov loop expectation values ($) and($) as a func- 
tion of T and /t. 



B. Quark- meson parameters 

The four parameters of the QM model, i.e. <?, A, v and 
c, are chosen such that chiral symmetry is spontaneously 
broken in the vacuum and the cr-field develops a finite 
expectation value (a) = f v , where f n = 93 MeV is set to 
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the pion decay constant. Due to the pseudoscalar char- 
acter of the pions the corresponding expectation values 
vanish, (7?) = 0. 

The Yukawa coupling constant g is fixed by the con- 
stituent quark mass in the vacuum g = m q j f w . Using 
the partially conserved axial vector current (PCAC) re- 
lation the explicit symmetry breaking parameter c is de- 
termined by c = tti^/tt, where is the pion mass. 
The quartic coupling constant A is given by the sigma 
mass m a via the relation A = (m 2 — m^)/(2f 2 ). Finally, 
the parameter v 2 is found by minimizing the potential 
in radial direction, yielding v 2 = (a) 2 - c/(A(<t». For 
the ground state where (a) = f n this expression can be 
rewritten as v 2 = f 2 —m 2 / A . It is positive in the Nambu- 
Goldstone phase. 

In the vacuum we fix the model parameters to = 



138 MeV, m a = 600 McV, U = 93 MeV and m„ 



300 



McV which result in c - 1.77 ■ 10 6 MeV 3 , v ~ 87.6 MeV, 
A - 19.7 and g ~ 3.2. 



C. Phase structure 

The phase structure of the PQM model is determined 
by the behavior of the order parameters a, $ and $ and 
of the grand canonical potential as a function of tem- 
perature and quark chemical potential. All numerical 
results have been obtained for Nf = 2. Then To = 208 
MeV in agreement with Tab. HI This value is different 
from that taken in Ref. 0, \bf$ where T = 270 McV, 
the value of Nf = 0. In these works T = 210 MeV 
has been fixed in order to compare with lattice results. 
The 7V/-depcndcncc suggested in the present work offers 
a qualitative explanation for this choice. 

In Fig.[T]the temperature dependence of the chiral con- 
densate (qq) and the Polyakov loop expectation value <E>, 
see Eq. (|2.4[) . at fjb = is shown in relative units. 



For vanishing chemical potential we have $ = <& as al- 
ready discussed. For T->oowe find $ ~ 1.11. Since the 
properly normalized expectation value <& tends towards 
unity we have normalized the mean fields accordingly. 
For temperatures at about 200-300 MeV the normalized 
3> has a /i-dependent maximum and decreases for larger 
temperatures towards one, see Fig. [2] This is in qualita- 
tive agreement with perturbation theory which predicts 
an increasing $ within an expansion around vanishing 
gauge fields, e.g. 



57 



58] 
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FIG. 1: The normalized chiral condensate (qq) and the 
Polyakov loop $ as a function of temperature for fi = 0. A 
chiral crossover is found at T ~ 180 MeV and a deconfinement 
crossover at a similar temperature. 




H = MeV 
H = 100 MeV 
ji = 168 MeV 
\l = 270 MeV 



At fi — we find a chiral crossover temperature 



200 400 600 800 1000 1200 1400 
T [MeV] 



FIG. 2: The normalized Polyakov loop variable "I> for large 
temperatures for several chemical potentials /1. 



T c = 184 MeV with an error of ~ ±14 MeV originat- 
ing in the error estimate ±30 MeV for To. For example, 
using the two-loop running of the coupling (|2.15p . and 
hence T (N f ) = 192 MeV we are led to T c - 177 MeV. 
In the presence of dynamical quarks the Polyakov loop 
shows also a crossover at the same pseudo-critical tem- 
perature. This can be read off from the peak position 
of d(qq)/dT and d$/dT. In Fig. [3] these quantities are 
shown as function of the temperature. 

In two-flavor lattice simulations extrapolated to the 
chiral limit a pseudo-critical temperature T c = 173 ± 8 
McV is found using improved staggered fermions [H . Re- 
cently, a recalculation of the transition temperature with 
staggered fermions for two light and one heavier quark 
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1.2 



120 140 160 180 200 220 240 
T [MeV] 

FIG. 3: The temperature dependence of d(qq)/dT and 9$/<9T 
for n = 0. The Polyakov variable is scaled by a factor of 5. 



mass close to their physical values yields a T c = 192 ± 7 
MeV using the Sommer parameter ro for the continuum 
extrapolation [5^] . This result has to be contrasted with 
another recent lattice analysis with staggered fermions 
but using four different sets of lattice sizes iV T = 4, 6, 8 
and 10 to perform the continuum extrapolation [601 ] . 
From the same physical observable this group finds a crit- 
ical temperature T c = 151 ±3 MeV. Functional RG stud- 

30, 



172±*° MeV 



ies yield a critical value of T c 
where the error originates in an estimate of the uncer- 
tainty similar to the considerations put forward here. On 
the other hand, using the same parameters for the quark- 
meson model without the Polyakov loop modifications a 
crossover temperature of T c ~ 150 MeV emerges [r3l| . 
This situation calls for refined studies both on the lat- 
tice as well as within functional methods to resolve the 
apparent quantitative inaccuracies. 

For finite /i the degeneracy of $ and $ disappears. 
The corresponding order parameters as function of tem- 
perature for several chemical potentials are collected in 
Figs. [4] and [5] For finite /i the Polyakov loop <t is always 
larger than It has a positive slope <9$/cfyt > for all 
temperatures, and peaks at some high temperature, see 
Fig. [21 Both, <I> and i> tend towards one for T — > 00. 

Above a critical chemical potential /j, c = 168 MeV all 
order parameters jump at the same temperature which 
signals a first-order phase transition. The critical end 
point (CEP) is found at (T c ,/i c ) = (150, 168) MeV. The 
corresponding chiral phase diagram obtained for a /1- 
independent To(Nf), (|2.17j) . is shown in Fig. [6] (upper 
lines). At the critical point the chiral first-order transi- 
tion line terminates and the transition becomes sccond- 



H=0 MeV 
H=168 MeV 
H=270 MeV 




150 200 
T [MeV] 



350 



FIG. 4: The normalized chiral quark condensate (<r) as a 
function of temperature for three different chemical potentials 
/i = 0, 168, 270 MeV. For /i = 270 MeV a first-order transition 
is found at T c ~ 81 MeV. 
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FIG. 5: Same as Fig. [4] for the normalized Polyakov loops $ 
and $. 



order, which induces a divergent quark number suscepti- 
bility. Lattice simulations are not conclusive concerning 
the existence and location of the critical point (ll [r32 . 67 1 . 

There are indications from lattice simulations at finite 
chemical potential that deconfinement and chiral symme- 
try restoration appear along the same critical line in the 
phase diagram. For the PQM model and //-independent 
To(Nf) the coincidence of deconfinement and chiral tran- 
sition at fx = disappears for finite fx. The deconfinement 
temperature is larger than the corresponding chiral tran- 
sition temperature. This is an unphysical scenario be- 
cause the deconfinement temperature should be smaller 
or equal to the chiral transition temperature. When 
resorting to the /i-dcpendcnt T (fi,Nf), (|2.22[) . we find 
coinciding transition lines for the entire phase diagram 
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50 100 150 200 250 300 350 
H [MeV] 

FIG. 6: Chiral phase diagrams for the PQM model. Upper 
lines for a /^-independent Polyakov loop potential and lower 
lines with /t-dependent corrections. The CEP's are approxi- 
mately located at (T c , fi c ) — (163, 164) MeV (upper case) and 
at (150, 168) MeV (lower case). 

within an accuracy of ±5 MeV. For this case the unique 
transition line lies slightly below the chiral one for the 
/i-independent choice To(Nf). This is shown in Fig. [6j 

D. Thermodynamic observables 

In order to investigate the influence of the Polyakov 
loop on the equilibrium thermodynamics we calculate 
several thermodynamic observables. All information of 
the system is contained in the grand canonical potential 
which is given by fJ in (|3 . 23[) evaluated at the mean-held 
level. 

We begin our analysis with the pressure of the system 
p. It is defined as the negative of the grand canonical 
potential and is normalized to vanish at T = fi = 0. 
In Fig. [7] the pressure divided by the QCD Stefan- 
Boltzmann (SB) limit is shown as function of the tem- 
perature for three different quark chemical potentials. 
The values of the chemical potentials are chosen such 
that one curve runs through the critical end point (CEP) 
(fi c = 168 MeV) and another curve through a first-order 
phase transition (/j, = 270 MeV). The QCD pressure in 
the SB limit for Nf massless quarks and (AT? — 1) mass- 
less gluons, relevant for the deconfined phase, is given 
by 



PSB 



{N 2 c -l)-+N c N f 



"7tt 2 


-I 




180" 







J_/£V 

12tt 2 \tJ 
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u = MeV 
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(3.29) 

where the first term denotes the gluonic contribution 



100 150 200 250 300 350 
T [MeV] 

FIG. 7: Scaled pressure p/pss for three different quark chem- 
ical potentials, = 0, 168, 270 MeV. T c {ji = 0) = 184 MeV. 



and the rest involves the fermions. The pressure is sup- 
pressed in the confined phase and starts to rise when 
deconfincment sets in. For all T and /i the pressure p/T 4 
stays below the QCD SB limit, a feature that is also ob- 
served in lattice calculations and other non-perturbativc 
approaches. For vanishing chemical potential the pres- 
sure is a smooth function of the temperature consistent 
with a crossover transition. At temperatures of twice the 
critical temperature the pressure reaches approximately 
80% of the SB limit. On the lattice two classes of data for 
the pressure obtained with a temporal extent N T = 4 and 
N T = 6 at /i = are currently available both of which are 
not extrapolated to the continuum 63j, [6_4| . Our results 
are in agreement with lattice simulations with a temporal 
extent of N T = 6 which is also closer to the continuum 
limit. This behavior is demonstrated in Fig.[5J 

An increase of the chemical potential leads to an in- 
crease of the pressure as more quark degrees of freedom 
are active. For a certain chemical potential the crossover 
transition changes to a first-order phase transition. In 
this case the pressure has a kink at the transition point 
but still is a continuous function. The kink at T ~ 100 
MeV for the fx = 270 MeV curve is clearly visible Fig. [7] 

At a first-order phase transition a finite latent heat 
builds up. This results in a jump of the entropy density 
s, which is defined as the negative derivative of the grand 
canonical potential with respect to the temperature. It 
is identical to the temperature derivative of the pressure. 
In Fig. [9] we show s divided by the corresponding QCD 
SB limit for the same chemical potentials as in the pre- 
ceding figure. The low- and high-temperature behavior 
of this quantity can be understood in a similar fashion as 
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FIG. 8: Scaled pressure p/psb for p = 0. The PQM model 
prediction (solid line) is compared to lattice results for N T — 4 
and N T — 6. Lattice data taken from Ref. [631 . 
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FIG. 10: Scaled pressure difference Ap/T 4 for three dif- 
ferent chemical potentials. The curves correspond to fi = 
100, 168, 270 MeV from below. 
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FIG. 9: Same as described in the legend to Fig. [7] for the 
scaled entropy s/ssb- 



those of the pressure. It is continuous in the vicinity of 
the crossover transition and reaches less than 40% of the 
SB limit around these temperatures. For chemical po- 
tential values larger than the critical one a finite latent 
heat emerges which further increases with the chemical 
potential. 

Another quantity that is accessible in lattice QCD at 
finite chemical potential is the pressure difference Ap. 
It is defined as Ap(T, /i) = p(T, /i) — p(T, /i = 0) and 
is Taylor-expanded around fi — in powers of the di- 
mensionless quantity fi/T on the lattice. Because odd 
derivatives of the free energy with respect to fi vanish 
only even powers appear in this expansion. In our model 
we have computed the pressure difference without refer- 
ring to an expansion. In Fig.[TO]the scaled pressure differ- 



ence Ap(T, /i)/T 4 versus temperature for three chemical 
potential values is shown. The bottom curve corresponds 
to /i = 100 MeV. It is always a continuous function and 
shows a kink at a first-order phase transition. Ap rises 
steeply across the chiral transition and peaks almost at 
the same temperature for all chemical potentials. For 
larger temperatures it decreases almost as 1/T 2 which 
follows from the SB limit. Nevertheless, due to the T- 
and //-dependent Polyakov fields slight deviations of the 
1/T 2 SB behavior are seen. 

Another interesting observable is the net quark den- 
sity. It is obtained from the thermodynamic potential 
via n q = —dfl(T,fi)/dfj.. The quark density, normal- 
ized to 1/T 3 , is displayed as a function of the tem- 
perature in Fig. [IT] for three different chemical poten- 
tials fi = 100, 168 and 270 MeV. In comparison to the 
pure quark-meson model without the Polyakov loop the 
quark density in the confined phase is much more sup- 
pressed when the interaction of quarks with the Polyakov 
loop is adde d I40L |61|. A similar effect is seen in the 
PNJL model [20(. Above the phase transition, the quark 
density of the pure quark-meson model approaches the 
Stefan-Boltzmann limit n q = Nf[i(T 2 + {^jl/tt) 2 ) imme- 
diately. With the Polyakov loop dynamics this behav- 
ior is changed drastically. The quark densities increase 
slightly above the corresponding SB limits and decrease 
again with growing temperature. For high temperatures 
the SB limit of the quark density is always reached from 
above. At a first-order phase transition n q jumps and 
drops immediately after the transition for increasing tem- 
peratures. 
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T [MeV] 

FIG. 11: Same as described in the legend to Fig. [10] for the 
quark number density n q /T 3 . The dashed lines denote the 
corresponding Stefan-Boltzmann limits. 

The quark number susceptibility measures the static 
response of the quark number density to an infinitesimal 
variation of the quark chemical potential and is given 
by Xq = driq/dfj,. It is shown in Fig. [T^] as a function 
of temperature for several /i. This observable can be 



H = MeV 
H= 168 MeV 
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FIG. 12: The scaled chiral susceptibility Xq/T 2 as a function 
of temperature for three different quark chemical potentials, 
fi = 0, 168, 270 MeV. 

used to verify the existence and location of the critical 
end point in the phase diagram. At a first-order phase 
transition this quantity has a discontinuity and only at a 
second-order critical end point it is divergent. Even for 
finite pion masses the critical point is of second-order and 
induces a divergent quark number susceptibility. This 
behavior is seen in Fig. [12] For fi = 168 MeV, close to 
the critical chemical potential, Xq diverges at the critical 
temperature. 



The modifications caused by the quark-gluon interac- 
tion on the quark number susceptibility, are similar as 
those already discussed in the context of the quark num- 
ber density. Compared to the pure quark-meson model 
Xq is again more suppressed below the chiral phase tran- 
sition. Above the transition \q ues above the correspond- 
ing SB limit Xq/T 2 = iV/(l + 3/7r 2 (^/T) 2 ). At high tem- 
peratures the SB limit (not shown in the figure) is again 
reached from above. 



IV. CONCLUSION 

In the present paper we have extended the Nf = 2 
quark-meson model to include certain aspects of gluon 
dynamics via the Polyakov loop. This PQM model com- 
bines important symmetry aspects in the limit of in- 
finitely heavy quarks with those of the light quark sec- 
tor of nearly massless up- and down quarks. Within the 
mean-field approximation we have discussed the ensuing 
phase diagram of strongly interacting matter in the grand 
canonical ensemble, as was done previously in a similar 
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65] . One of the ben- 



extension of the NJL model |39|, 
efits is an improvement of the thermodynamical behavior 
of several bulk quantities such as pressure, entropy etc. 
when compared to lattice data. As a novelty we propose 
to include the Nf and \i dependence of the the running 
coupling a in the parameter of the Polyakov loop poten- 
tial. A qualitative estimate is provided by the one loop 
/3-function for the gauge coupling a as well as using the 
hard dense loop approximation. Then we are led to a 
Nf and fi dependent To, the critical temperature of the 
Polyakov loop model, which decreases with increasing Nf 
and fi. These modifications already involve coinciding 
peaks in the temperature derivative of the Polyakov loop 
expectation value and the chiral condensate at fx = 0, 
in agreement with the lattice findings of 66, 6^]. Inter- 
estingly this coincidence of the deconfinement and chiral 
symmetry restoration persists at finite [i. 

The findings of the present work provide a promising 
startin g p oint for a functional RG study in the present 
model [26(, and further extensions towards full QCD: in 
particular, we aim at removing the perturbative nature 
of the above estimates as well as allowing for a fully cou- 
pled PQM model. The last step consists of including the 
full gauge dynamics. This is particularly relevant for the 
important issue of the Polyakov loop potential at finite 

being intimately related to the open question of the 
existence and location of the critical point in the QCD 
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phase diagram. 
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